Analysis of Immune Intratumor Heterogeneity Highlights Immunoregulatory and Coinhibitory Lymphocytes as Hallmarks of Recurrence in Stage I Non-Small Cell Lung Cancer.

Our understanding of the molecular mechanisms underlying postsurgical recurrence of non-small cell lung cancer (NSCLC) is rudimentary. Molecular and T cell repertoire intratumor heterogeneity (ITH) have been reported to be associated with postsurgical relapse; however, how ITH at the cellular level impacts survival is largely unknown. Here we report the analysis of 2880 multispectral images representing 14.2% to 27% of tumor areas from 33 patients with stage I NSCLC, including 17 cases (relapsed within 3 years after surgery) and 16 controls (without recurrence ≥5 years after surgery) using multiplex immunofluorescence. Spatial analysis was conducted to quantify the minimum distance between different cell types and immune cell infiltration around malignant cells. Immune ITH was defined as the variance of immune cells from 3 intratumor regions. We found that tumors from patients having relapsed display different immune biology compared with nonrecurrent tumors, with a higher percentage of tumor cells and macrophages expressing PD-L1 (P =.031 and P =.024, respectively), along with an increase in regulatory T cells (Treg) (P =.018), antigen-experienced T cells (P =.025), and effector-memory T cells (P =.041). Spatial analysis revealed that a higher level of infiltration of PD-L1+ macrophages (CD68+PD-L1+) or antigen-experienced cytotoxic T cells (CD3+CD8+PD-1+) in the tumor was associated with poor overall survival (P =.021 and P =.006, respectively). A higher degree of Treg ITH was associated with inferior recurrence-free survival regardless of tumor mutational burden (P =.022), neoantigen burden (P =.021), genomic ITH (P =.012) and T cell repertoire ITH (P =.001). Using multiregion multiplex immunofluorescence, we characterized ITH at the immune cell level along with whole exome and T cell repertoire sequencing from the same tumor regions. This approach highlights the role of immunoregulatory and coinhibitory signals as well as their spatial distribution and ITH that define the hallmarks of tumor relapse of stage I NSCLC.

Immunogenomic intertumor heterogeneity across primary and metastatic sites in a patient with lung adenocarcinoma.

BACKGROUND: Lung cancer is the leading cause of cancer death, partially owing to its extensive heterogeneity. The analysis of intertumor heterogeneity has been limited by an inability to concurrently obtain tissue from synchronous metastases unaltered by multiple prior lines of therapy. METHODS: In order to study the relationship between genomic, epigenomic and T cell repertoire heterogeneity in a rare autopsy case from a 32-year-old female never-smoker with left lung primary late-stage lung adenocarcinoma (LUAD), we did whole-exome sequencing (WES), DNA methylation and T cell receptor (TCR) sequencing to characterize the immunogenomic landscape of one primary and 19 synchronous metastatic tumors. RESULTS: We observed heterogeneous mutation, methylation, and T cell patterns across distinct metastases. Only TP53 mutation was detected in all tumors suggesting an early event while other cancer gene mutations were later events which may have followed subclonal diversification. A set of prevalent T cell clonotypes were completely excluded from left-side thoracic tumors indicating distinct T cell repertoire profiles between left-side and non left-side thoracic tumors. Though a limited number of predicted neoantigens were shared, these were associated with homology of the T cell repertoire across metastases. Lastly, ratio of methylated neoantigen coding mutations was negatively associated with T-cell density, richness and clonality, suggesting neoantigen methylation may partially drive immunosuppression. CONCLUSIONS: Our study demonstrates heterogeneous genomic and T cell profiles across synchronous metastases and how restriction of unique T cell clonotypes within an individual may differentially shape the genomic and epigenomic landscapes of synchronous lung metastases.

Expression of the SAR2-Cov-2 receptor ACE2 reveals the susceptibility of COVID-19 in non-small cell lung cancer.

Recent studies have revealed that cancer patients had a higher risk of having coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), compared to patients without cancer. The expression of angiotensin-converting enzyme 2 (ACE2), the receptor of SARS-CoV-2, was aberrantly expressed in many tumors. In this study, by exploring the TCGA and GTEx public databases, we investigated ACE2 expression and its association with prognosis in non-small cell lung cancer (NSCLC), the most susceptible caner type. We found that lung was one of the major organs with highly expressed ACE2. Furthermore, ACE2 expression was significantly elevated in lung adenocarcioma (LUAD) and lung squamous cell carcinoma (LUSC) compared to normal tissues. DNA methylation might be one possible mechanism leading to ACE2 upregulation. Despite that, the AEC2 expression was not statistically associated with disease-free survival (DFS) and overall survival (OS) for LUAD patients, and higher ACE2 expression was associated with prolonged DFS in LUSC patients. Taken together, we observed ACE2 was highly expressed in LUAD and LUSC despite the controversial role of ACE2 expression in predicting prognosis in these two common lung cancer types.

Integrated Bioinformatics Analysis Reveals Key Candidate Genes and Pathways Associated With Clinical Outcome in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) accounts for approximately 85-90% of all liver cancer cases and has poor relapse-free survival. There are many gene expression studies that have been performed to elucidate the genetic landscape and driver pathways leading to HCC. However, existing studies have been limited by the sample size and thus the pathogenesis of HCC is still unclear. In this study, we performed an integrated characterization using four independent datasets including 320 HCC samples and 270 normal liver tissues to identify the candidate genes and pathways in the progression of HCC. A total of 89 consistent differentially expression genes (DEGs) were identified. Gene-set enrichment analysis revealed that these genes were significantly enriched for cellular response to zinc ion in biological process group, collagen trimer in the cellular component group, extracellular matrix (ECM) structural constituent conferring tensile strength in the molecular function group, protein digestion and absorption, mineral absorption and ECM-receptor interaction. Network system biology based on the protein-protein interaction (PPI) network was also performed to identify the most connected and important genes based on our DEGs. The top five hub genes including osteopontin (SPP1), Collagen alpha-2(I) chain (COL1A2), Insulin-like growth factor I (IGF1), lipoprotein A (LPA), and Galectin-3 (LGALS3) were identified. Western blot and immunohistochemistry analysis were employed to verify the differential protein expression of hub genes in HCC patients. More importantly, we identified that these five hub genes were significantly associated with poor disease-free survival and overall survival. In summary, we have identified a potential clinical significance of these genes as prognostic biomarkers for HCC patients who would benefit from experimental approaches to obtain optimal outcome.

Genomic assessment distinguishes intrapulmonary metastases from synchronous primary lung cancers.

BACKGROUND: Multiple synchronous lung tumors (MSLT), particularly within a single lobe, represent a diagnostic and treatment challenge. While histologic assessment was once the only method to possibly distinguish multiple primary lung cancers, there is a growing interest in identifying unique genomic features or mutations to best characterize these processes. METHODS: In order to differentiate multiple primary lung malignancies from intrapulmonary metastases in patients with MSLT, we performed whole exome sequencing (WES) on 10 tumor samples from 4 patients with MSLT. RESULTS: Shared mutations between tumors from the same patient varied from 0-91%. Patient 3 shared no common mutations; however, in Patients 2 and 4, identical mutations were identified among all tumors from each patient, suggesting that the three tumors identified in Patient 3 represent separate primary lung cancers, while those of Patients 1, 2 and 4 signify hematogenous and lymphatic spread. CONCLUSIONS: A high proportion of shared mutations between different lung tumors is likely indicative of intrapulmonary metastatic disease, while tumors with distinct genomic profiles likely represent multiple primary malignancies driven by distinct molecular events. Application of genomic profiling in the clinical setting may prove to be important to precise management of patients with MSLT.

Comprehensive T cell repertoire characterization of non-small cell lung cancer.

Immunotherapy targeting T cells is increasingly utilized to treat solid tumors including non-small cell lung cancer (NSCLC). This requires a better understanding of the T cells in the lungs of patients with NSCLC. Here, we report T cell repertoire analysis in a cohort of 236 early-stage NSCLC patients. T cell repertoire attributes are associated with clinicopathologic features, mutational and immune landscape. A considerable proportion of the most prevalent T cells in tumors are also prevalent in the uninvolved tumor-adjacent lungs and appear specific to shared background mutations or viral infections. Patients with higher T cell repertoire homology between the tumor and uninvolved tumor-adjacent lung, suggesting a less tumor-focused T cell response, exhibit inferior survival. These findings indicate that a concise understanding of antigens and T cells in NSCLC is needed to improve therapeutic efficacy and reduce toxicity with immunotherapy, particularly adoptive T cell therapy.

Correction: Distinct co-acquired alterations and genomic evolution during TKI treatment in non-small-cell lung cancer patients with or without acquired T790M mutation.

The original version of this Article omitted the following from the Acknowledgements: Professor Stebbing sits on SABs for Celltrion, Singapore Biotech, Vor Biopharma, TLC Biopharmaceuticals and Benevolent AI, has consulted with Lansdowne partners, Vitruvian and Social Impact Capital and Chairs the Board of Directors for BB Biotech Healthcare Trust and Xerion Healthcare. This has now been corrected in both the PDF and HTML versions of the Article.

Distinct co-acquired alterations and genomic evolution during TKI treatment in non-small-cell lung cancer patients with or without acquired T790M mutation.

EGFR-mutant non-small-cell lung cancer (NSCLC) patients inevitably develop drug resistance when treated with EGFR tyrosine kinase inhibitors (TKIs). Systematic genetic analysis is important to understand drug-resistant mechanisms; however, the clinical significance of co-occurring genetic alterations at baseline, co-acquired mutations at progressive disease (PD), and the clonal evolution remain underinvestigated. We performed targeted sequencing of pre-treatment and PD tumor samples from 54 EGFR-mutant NSCLC patients. Ten additional patients were sequenced using whole-exome sequencing to infer the clonal evolution patterns. We observed a domain-dependent effect of PIK3CA mutation at baseline on patient progression-free survival (PFS). In addition, at baseline, 9q34.3/19p13.3 (NOTCH1/STK11/GNA11) showed a co-deletion pattern, which was associated with a significantly worse PFS (p = 0.00079). T790M-postive patients with other concurrent acquired oncogenic mutations had a significantly shorter PFS (p = 0.005). Besides acquired T790M mutation, chromosomal instability (CIN) related genes, including AURKA and TP53 alterations, were the most frequently acquired events. CIN significantly increased during TKI treatment in T790M-negative patients and is a candidate resistance mechanism to the first-generation TKIs. Clonal evolution analyses suggest that the composition and relationship among resistant subclones, particularly relationship with T790M subclone, affect patients' outcomes. Overall, our findings of novel co-occurring alterations and clonal evolution patterns can be served as predictive biomarkers to stratify patients and help to better understand the drug-resistant mechanism to TKIs.

Co-fuse: a new class discovery analysis tool to identify and prioritize recurrent fusion genes from RNA-sequencing data.

Recurrent oncogenic fusion genes play a critical role in the development of various cancers and diseases and provide, in some cases, excellent therapeutic targets. To date, analysis tools that can identify and compare recurrent fusion genes across multiple samples have not been available to researchers. To address this deficiency, we developed Co-occurrence Fusion (Co-fuse), a new and easy to use software tool that enables biologists to merge RNA-seq information, allowing them to identify recurrent fusion genes, without the need for exhaustive data processing. Notably, Co-fuse is based on pattern mining and statistical analysis which enables the identification of hidden patterns of recurrent fusion genes. In this report, we show that Co-fuse can be used to identify 2 distinct groups within a set of 49 leukemic cell lines based on their recurrent fusion genes: a multiple myeloma (MM) samples-enriched cluster and an acute myeloid leukemia (AML) samples-enriched cluster. Our experimental results further demonstrate that Co-fuse can identify known driver fusion genes (e.g., IGH-MYC, IGH-WHSC1) in MM, when compared to AML samples, indicating the potential of Co-fuse to aid the discovery of yet unknown driver fusion genes through cohort comparisons. Additionally, using a 272 primary glioma sample RNA-seq dataset, Co-fuse was able to validate recurrent fusion genes, further demonstrating the power of this analysis tool to identify recurrent fusion genes. Taken together, Co-fuse is a powerful new analysis tool that can be readily applied to large RNA-seq datasets, and may lead to the discovery of new disease subgroups and potentially new driver genes, for which, targeted therapies could be developed. The Co-fuse R source code is publicly available at https://github.com/sakrapee/co-fuse .

TCR Repertoire Intratumor Heterogeneity in Localized Lung Adenocarcinomas: An Association with Predicted Neoantigen Heterogeneity and Postsurgical Recurrence.

Genomic intratumor heterogeneity (ITH) may be associated with postsurgical relapse of localized lung adenocarcinomas. Recently, mutations, through generation of neoantigens, were shown to alter tumor immunogenicity through T-cell responses. Here, we performed sequencing of the T-cell receptor (TCR) in 45 tumor regions from 11 localized lung adenocarcinomas and observed substantial intratumor differences in T-cell density and clonality with the majority of T-cell clones restricted to individual tumor regions. TCR ITH positively correlated with predicted neoantigen ITH, suggesting that spatial differences in the T-cell repertoire may be driven by distinct neoantigens in different tumor regions. Finally, a higher degree of TCR ITH was associated with an increased risk of postsurgical relapse and shorter disease-free survival, suggesting a potential clinical significance of T-cell repertoire heterogeneity.Significance: The present study provides insights into the ITH of the T-cell repertoire in localized lung adenocarcinomas and its potential biological and clinical impact. The results suggest that T-cell repertoire ITH may be tightly associated to genomic ITH and disease relapse. Cancer Discov; 7(10); 1088-97. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 1047.

DNA methylation intratumor heterogeneity in localized lung adenocarcinomas.

Cancers are composed of cells with distinct molecular and phenotypic features within a given tumor, a phenomenon termed intratumor heterogeneity (ITH). Previously, we have demonstrated genomic ITH in localized lung adenocarcinomas; however, the nature of methylation ITH in lung cancers has not been well investigated. In this study, we generated methylation profiles of 48 spatially separated tumor regions from 11 localized lung adenocarcinomas and their matched normal lung tissues using Illumina Infinium Human Methylation 450K BeadChip array. We observed methylation ITH within the same tumors, but to a much less extent compared to inter-individual heterogeneity. On average, 25% of all differentially methylated probes compared to matched normal lung tissues were shared by all regions from the same tumors. This is in contrast to somatic mutations, of which approximately 77% were shared events amongst all regions of individual tumors, suggesting that while the majority of somatic mutations were early clonal events, the tumor-specific DNA methylation might be associated with later branched evolution of these 11 tumors. Furthermore, our data showed that a higher extent of DNA methylation ITH was associated with larger tumor size (average Euclidean distance of 35.64 (> 3cm, median size) versus 27.24 (<= 3cm), p = 0.014), advanced age (average Euclidean distance of 34.95 (above 65) verse 28.06 (below 65), p = 0.046) and increased risk of postsurgical recurrence (average Euclidean distance of 35.65 (relapsed patients) versus 29.03 (patients without relapsed), p = 0.039).

Antimicrobial susceptibility and clarithromycin resistance patterns of Helicobacter pylori clinical isolates in Vietnam.

Helicobacter pylori is a gastric pathogen that causes several gastroduodenal disorders such as peptic ulcer disease and gastric cancer.  Eradication efforts of H. pylori are often hampered by antimicrobial resistance in many countries, including Vietnam.  Here, the study aimed to investigate the occurrence of antimicrobial resistance among H. pylori clinical isolates across 13 hospitals in Vietnam.  The study further evaluated the clarithromycin resistance patterns of H. pylori strains.  In order to address the study interests, antimicrobial susceptibility testing, epsilometer test and PCR-based sequencing were performed on a total of 193 strains isolated from patients, including 136 children (3-15 years of age) and 57 adults (19-69 years of age).  Antimicrobial susceptibility testing showed that the overall resistance to amoxicillin, clarithromycin, levofloxacin, metronidazole, and tetracycline was 10.4%, 85.5%, 24.4%, 37.8%, and 23.8% respectively.  The distribution of minimum inhibitory concentrations (MICs) of clarithromycin-resistant strains was 85.5% with MIC >0.5 µg/mL.  The majority of the clarithromycin resistant isolates (135 of 165 subjects) have MICs ranging from 2 µg/mL to 16 µg/mL.  Furthermore, sequencing detection of mutations in 23S rRNA gene revealed that strains resistant and susceptible to clarithromycin contained both A2143G and T2182C mutations.  Of all isolates, eight clarithromycin-resistant isolates (MIC >0.5 µg/mL) had no mutations in the 23S rRNA gene.  Collectively, these results demonstrated that a proportion of clarithromycin-resistant H. pylori strains, which are not related to the 23S rRNA gene mutations, could be potentially related to other mechanisms such as the presence of an efflux pump or polymorphisms in the CYP2C19 gene.  Therefore, the present study suggests that providing susceptibility testing prior to treatment or alternative screening strategies for antimicrobial resistance is important for future clinical practice.  Further studies on clinical guidelines and treatment efficacy are pivotal for successful eradication of H. pylori infection.

Genomic analyses identify molecular subtypes of pancreatic cancer.

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-ß, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Whole-genome characterization of chemoresistant ovarian cancer.

Patients with high-grade serous ovarian cancer (HGSC) have experienced little improvement in overall survival, and standard treatment has not advanced beyond platinum-based combination chemotherapy, during the past 30 years. To understand the drivers of clinical phenotypes better, here we use whole-genome sequencing of tumour and germline DNA samples from 92 patients with primary refractory, resistant, sensitive and matched acquired resistant disease. We show that gene breakage commonly inactivates the tumour suppressors RB1, NF1, RAD51B and PTEN in HGSC, and contributes to acquired chemotherapy resistance. CCNE1 amplification was common in primary resistant and refractory disease. We observed several molecular events associated with acquired resistance, including multiple independent reversions of germline BRCA1 or BRCA2 mutations in individual patients, loss of BRCA1 promoter methylation, an alteration in molecular subtype, and recurrent promoter fusion associated with overexpression of the drug efflux pump MDR1.

Whole genomes redefine the mutational landscape of pancreatic cancer.

Pancreatic cancer remains one of the most lethal of malignancies and a major health burden. We performed whole-genome sequencing and copy number variation (CNV) analysis of 100 pancreatic ductal adenocarcinomas (PDACs). Chromosomal rearrangements leading to gene disruption were prevalent, affecting genes known to be important in pancreatic cancer (TP53, SMAD4, CDKN2A, ARID1A and ROBO2) and new candidate drivers of pancreatic carcinogenesis (KDM6A and PREX2). Patterns of structural variation (variation in chromosomal structure) classified PDACs into 4 subtypes with potential clinical utility: the subtypes were termed stable, locally rearranged, scattered and unstable. A significant proportion harboured focal amplifications, many of which contained druggable oncogenes (ERBB2, MET, FGFR1, CDK6, PIK3R3 and PIK3CA), but at low individual patient prevalence. Genomic instability co-segregated with inactivation of DNA maintenance genes (BRCA1, BRCA2 or PALB2) and a mutational signature of DNA damage repair deficiency. Of 8 patients who received platinum therapy, 4 of 5 individuals with these measures of defective DNA maintenance responded.

Genomic catastrophes frequently arise in esophageal adenocarcinoma and drive tumorigenesis.

Oesophageal adenocarcinoma (EAC) incidence is rapidly increasing in Western countries. A better understanding of EAC underpins efforts to improve early detection and treatment outcomes. While large EAC exome sequencing efforts to date have found recurrent loss-of-function mutations, oncogenic driving events have been underrepresented. Here we use a combination of whole-genome sequencing (WGS) and single-nucleotide polymorphism-array profiling to show that genomic catastrophes are frequent in EAC, with almost a third (32%, n=40/123) undergoing chromothriptic events. WGS of 22 EAC cases show that catastrophes may lead to oncogene amplification through chromothripsis-derived double-minute chromosome formation (MYC and MDM2) or breakage-fusion-bridge (KRAS, MDM2 and RFC3). Telomere shortening is more prominent in EACs bearing localized complex rearrangements. Mutational signature analysis also confirms that extreme genomic instability in EAC can be driven by somatic BRCA2 mutations. These findings suggest that genomic catastrophes have a significant role in the malignant transformation of EAC.

A workflow to increase verification rate of chromosomal structural rearrangements using high-throughput next-generation sequencing.

Somatic rearrangements, which are commonly found in human cancer genomes, contribute to the progression and maintenance of cancers. Conventionally, the verification of somatic rearrangements comprises many manual steps and Sanger sequencing. This is labor intensive when verifying a large number of rearrangements in a large cohort. To increase the verification throughput, we devised a high-throughput workflow that utilizes benchtop next-generation sequencing and in-house bioinformatics tools to link the laboratory processes. In the proposed workflow, primers are automatically designed. PCR and an optional gel electrophoresis step to confirm the somatic nature of the rearrangements are performed. PCR products of somatic events are pooled for Ion Torrent PGM and/or Illumina MiSeq sequencing, the resulting sequence reads are assembled into consensus contigs by a consensus assembler, and an automated BLAT is used to resolve the breakpoints to base level. We compared sequences and breakpoints of verified somatic rearrangements between the conventional and high-throughput workflow. The results showed that next-generation sequencing methods are comparable to conventional Sanger sequencing. The identified breakpoints obtained from next-generation sequencing methods were highly accurate and reproducible. Furthermore, the proposed workflow allows hundreds of events to be processed in a shorter time frame compared with the conventional workflow.